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ABSTRACT 

A qualitative study of the surface trapping of low-frequency non-radial g-modes in 
rotating early-type stars is undertaken within the Cowling, adiabatic and traditional 
approximations. A dispersion relation describing the local character of waves in a ro- 
tating star is derived; this dispersion relation is then used to construct propagation 
diagrams for a 7 Mq stellar model, which show the location and extent of wave trap- 
ping zones inside the star. It is demonstrated that, at frequencies below a cut-off, 
waves cannot be fully trapped within the star, and will leak through the surface. Ex- 
pressions for the cut-off frequency are derived in both the non-rotating and rotating 
cases; it is found from these expressions that the cut-off frequency increases with the 
rotation rate for all but prograde sectoral modes. 

While waves below the cut-off cannot be reflected at the stellar surface, the pres- 
ence of a sub-surface convective region in the stellar model, due to Hen ionization, 
means that they can become partially trapped within the star. The energy leakage 
associated with such waves, which are assigned the moniker virtual modes due to their 
discrete eigenfrequencies, means that stability analyses which disregard their existence 
(by assuming perfect reflection at the stellar surface) may be in error. 

The results are of possible relevance to the 53 Per and SPB classes of variable 
star, which exhibit pulsation frequencies of the same order of magnitude as the cut- 
off frequencies found for the stellar model. It is suggested that observations either of 
an upper limit on variability periods (corresponding to the cut-off), or of line-profile 
variations due to virtual modes, may permit asteroseismological studies of the outer 
layers of these systems. 
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1 INTRODUCTION 

The self-excitation of global non-radial pulsation modes in a 
star is a prime example of positive feedback, whereby small 
oscillatory perturbations grow in amplitude via the efficient 
conversion of heat into vibrational energy by a suitable driv- 



ing mechanism (see, e.g., Unno el al. 1989 for a eoinprelien- 
sive re view of the topic). A fundamental ingredient in the 

feedback loop is that the oscillations must be trapped in 
some part of the stellar interior, so that energy does not 
leak from the system faster than it can be generated. Such 
trapping can occur when a pair of evanescent regions, where 
traveling waves cannot be supported, enclose a propagative 
region; waves are repeatedly reflected at the two evanescent 
boundaries, and the resulting superposition leads to a stand- 
ing wave of the normal-mode type. 

For waves excited in stellar envelopes, it is common for 
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the surface layers to serve as one of the evanescent regions 
requ ired f or the formation of a trapping zone. Ando & Os- 
aki (1975) demonstrated that such a situation occurs in the 
Sun, where low-order p-modes are trapped beneath the pho- 
tosphere, supporting a model first put forward by Ulrich 
([L970|) to e xplain the five-minute solar oscillation (Lcighton 
ct al. 196$ . However, the trapping is only effective for modes 



with frequencies below some cut-off; higher-frequency modes 
cannot be reflected at the photosphere, and will leak through 
the stellar surfa ce. T his issue was addressed in detail by 
Ando & Osaki (1977), who found that, although leakage 
does occur through the solar photosphere at frequencies 
above the cut-off, some waves can subsequently be reflected 
at the chromosphere-corona interface, and standing wave s 
are able to form. More recently, Balmforth & Gough (199C) 



suggested that such coronal reflection can explain appar- 
ent obs ervations of high-freque ncy chromospheric standing 
waves (Fleck & Deubner 1989), although debate concern- 
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ing this interpretation still continues ([Kumar et al. 1994 



Dzhalilov fc Staude 1995; Jefferies 199^) 



Pulsation in massive, early stars (types O and B) is 
qualitatively quite different from the solar case, due to 
the gross structural differences between the two stellar 
classes. However, it is still subject to the same wave trap- 
ping requirements, since th e und erlying physics remains the 
same. Shibahashi & Osaki ( 197£| ), in their study of g-modes 
trapped within the hydrogen-burning shell of evolved mas- 
sive stars, found that high-frequency (low-order) modes can 
tunnel through an evanescent region separating core and en- 
velope, and thence escape from t he sta r. A complementary 
situation was discussed by Osaki (1977) when studying pul- 



sation in Cepheid-type stars; non-radial p-modes trapped 
within the envelope were able to tunnel through an evanes- 
cent region into the core, where they were damped rapidly 
without reflection at the centre. 

In both cases, the appropriate region of the star was 
modeled as an isolating oscillating unit with the inclusion 
of wave leakage at one boundary. The leakage was found 
to stabilize some modes which would otherwise have been 
self-excited, due to the asso ciated loss of vibrational energy 
from the star. Shibahashi (1979) analyzed wave trapping 



in an idealized stellar model (corresponding to an evolved 
massive star) using an asymptotic method, and discussed in 
some depth these two cases; in addition, he considered the 
situation where low-frequency (high-order) gr-modes are able 
to tunnel through an evanescent region in the envelope and 
thence escape through the stellar surface. 

More recently, however, relatively little attention has 
been shown regarding wave trapping issues at the surface 
of early- type stars; in particu l ar, stability analyses ( Cox et 
al. 1992J; [Kiriakidis et al. 199 2 ; Dziembowski & Pamyatnykh 



1995 



1 



aAitscliy*^^aio~1993| ) , based on the new opaci ty ca l- 



culations of Rogers & Iglesias (1992) and Seaton (1992) 



have assumed that the Lagrangian pressure perturbation Sp 
tends to zero or some limiting value at the stellar surface. 
Such an assumption corresponds to the ab initio condition 
that waves incident from the interior are totally reflected 
at the stellar surface; the possibility of leakage is thereby 
disregarded, and no consideration of trapping issues is un- 
dertaken. 

This is the first in a short series of papers studying 
the surface trapping of low-frequency p-modes in early-type 
stars, in an attempt to re-open discussion of, and investiga- 
tion into, this important area. Much of the w ork is conceptu- 
ally developed from that of Ando & Osaki ( 1975| ); however, 
in light of recent rese arch into the influence of rotation on 



low-frequency mo des (Lee & Saio 1990 



Bildsten et al. 1996 



Lee fc Saio 1997), and due to the fact that significant rota- 



tion a ppears to be commo nplace in O- and B-star popula- 
tions (Howarth et al. 1997), the theory is updated to include 
rotational effects. 

The current paper serves as a introduction, covering 
the more qualitative, general aspects of the study; subse- 
quent papers will investigate various issues arising from this 
paper in greater depth. The following section reviews the 
pulsation equations appropriate for low-frequency g-modes 
in rotating stars, whilst section [] derives the dispersion rela- 
tion corresponding to these equations. The trapping of waves 
described by this dispersion relation is examined in section 
W with the aid of propagation diagrams, and the effect of 



rotation on the eigenfrequencies of individual modes is dis- 
cussed in section g. The findings are discussed in section ^, 
and summarized in section [jj. 



2 PULSATION EQUATIONS 

The dynamics of pulsation in a rotating star differ from 
the non-rotating case due to the influence of the fictitious 
Coriolis and centrifugal forces, which arise as a consequence 
of the non-inertial nature of a rotating frame of reference. 
The centrifugal force breaks the equilibrium symmetry of 
the star, so that the level (equipotential) surfaces become 
oblate spheroids rather than the usual concentric spheres. 
Such a change in stellar configuration will manifest itself 
implicitly in the pulsation equations, through modifications 
to the equilibrium variables of state. However, in the case of 
uniform (solid-body) rotation, no explicit modification of the 
pulsation equation s occurs due to this centrifugal distortion 
( [Unno et al. 1989| ). In contrast, the Coriolis force enters the 
pulsation equations explicitly, through the introduction of a 
velocity-dependent term in the hydrodynamical momentum 
equation. This term can lead to the significant modification 
of individual pulsation modes, and is also respo nsible for 



the existence of new classes of wave-like solutions ( Longuet- 



Higgins 196S ) which are not found in non-rotating systems. 

Simultaneous treatment of both forces within a pulsa- 
tion framework is fraught with difficulty. Some progress to- 
wards this goal has been made (e.g., Lee 1993; Lee & Baraffc 
1995), but attempts remain frustrated by the fact that the 
centrifugal distortion cannot really be considered as an a 
posteriori modification to the structure of a given star, but 
must be treated self-consistently with the evolution of the 
star (see, e.g., Meynet & Maeder 1997). However, in a num- 
ber of limiting cases, certain approximations can be made 
which simplify the problem significantly. In the case of the 
low-frequency modes, the Coriolis force will dominate the 
centrifugal force, and the effects of the latter on the equi- 
librium configuration may be disregarded if the rotation is 
not too sever e. Furthermo re, the so-called 'traditional ap- 
proximation' ( Eckart 1960| ) may be employed, whereby the 
horizontal component of the angular frequency vector of ro- 
tation CI is neglected. This approximation is most appropri- 
ate for lo w-frequency puls ation modes in the outer regions 
of a star (Unno et al. 1989), and therefore can be considered 



useful in the present study. 

In combination with the Cowling ( 1941 ) and adiabatic 
approximations, where the perturbations to the gravita- 
tional potential and specific entropy, respectively, are ne- 
glected, the traditional approximation renders the pulsa- 
tion equations separable in the spherical polar co-ordinates 
(r, 9, (f>). Solutions for the dependent variables £ r and p', the 
radial fluid displacement and Eulerian pressure per turba- 
tion, resp ectively, may then be written in the form (Lee & 
Saio 1997j) 



£r = Cr(r)6r(M; v) exp[i(m<?!> + cot)], 
p = p'(r)er(/-i; v) exp[i(m0 + ut)], 



(1) 
(2) 



where fj, = cos 9 is the normalized latitudinal distance from 
the equatorial plane, to is the pulsation frequency in the 
co-ro tating reference fram e , and OY l (p,\v) is a Hough func- 
tion flBildsten et al. 1996| |Lee fc Saio 1997|) . These Hough 
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f unctions are the eigen solutions of Laplace's tidal equation 
( Longuct-Higgins 1968), and form a one-parameter family 
in v = 2Q/u>, where fl = |fi| is the angular frequency of ro- 
tation. The integer indices I and m, with I ^ and \m\ ^ I, 
correspond to the harmonic degree and azimuthal order, re- 
s pectively, of the associated Legendre polynomials Pj m (/i) 
(Abramowitz & Stegun 1964) to which the Hough functions 
reduce in the non-rotating limit, so that OJ n (/i; 0) = P™(p). 
This index ing scheme, based on the one adopted by Le e & 
Saio (199C), is less general than that of Lee & Saio ( 1997 ), in 
that it does not encompass the Hough functions correspond- 
ing to Rossby and oscillatory convective modes (which do 
not have non-rotating counterparts); however, such modes 
are not considered herein, and the current scheme is suffi- 
cient. Note that ui is considered to be positive throughout 
the following discussion, and, therefore, prograde and retro- 
grade modes correspond to negative and positive values of 
m, respectively. 

The radial dependence of the solutions (|l|-|^) is de- 
scribed by the eigenfunctions £ r (r) and p'(r), which are gov- 
erned by a pair of coupled first-order differential equations. 
In order to facilitate subsequent manipulation, it is useful 
to write these equations in the form 



2 \ P 
Ul — 



(3) 



(4) 



i d a 2 t \ g c i 

r 2 dr y ' c£ oj 2 < 
and 

1 dp' g p' , 2 M 2,, 
-— + — — = (w -N )£r, 
p dr ci p 

where p, c s , g and are the local equilibrium values of 
the density, adiabatic sound speed, gravitational accelera- 
tion and Brunt- Vaisala frequency, respectively. Note that £ r 
and p' are now taken to be functions of r alone in both these 
and subsequent equations, unless explicitly stated. 

The quantity A; m appearing in equation (H), which 
arises as separation constant when solutions of the form 
(^-^) are sought, is the eigenvalue of Laplace's tidal equation 
corresponding to the appropriate Hough function OJ™(/x; v)- 
In the limit v = 0, this eigenvalue is equal to /(/ + 1), and 
equations (^-^) are then identical to those appropriate for a 
non-rotating star (e.g., Unno et al. 1989, §15.1). The utility 
of the traditional approximation thus lies in the fact that 
much of the formalism of the non-rotating case may also 
be applied to rotating stars with the simple replac ement of 
1(1 + 1) by A; m , a result first found by Lee & Saio (1987a). 

Global solution of equations W) must typically be ap- 
proached numerically; however, an examination of the local 
character of the solutions suffices in the present qualitative 
context. This character is governed by the dispersion rela- 
tion applicable to the equations, discussed in the following 
section. 



3 DISPERSION RELATION 

To derive a local dispersion relation for the pulsation equa- 
tions (j^jj) , it is useful first to place the equations in a canon- 
ical form similar to that introduced by Osaki (1975) for the 
non-rotating case. By defining the two new eigenfunctions, 



£ = r £ r exp 



(5) 



>l - 



P 



■ exp 



N 2 



dr 



(6) 



the left-hand sides of both pulsation equations may be writ- 
ten as a single derivative, and the canonical form is found 
as 



dr ciLO 1 



dr\ 
dr 



A; m C s 



UJ 2 | T], 



2 h(r) 



N 2 ) I 



where 



h(r) 



exp 



N 2 
9 



dr 



(7) 



(8) 



0) 



Note that h(r) is always positive, so that the original eigen- 
functions £ r and p' everywhere share the same sign as £ and 
f), respectively. 

Qualitative solution of these canonical e quati ons is ac- 
complished using the same method as Osaki (1975), namely, 



by assuming that the coefficients on the right-hand sides 
are independent of r. Such an assumption will be valid if 
the characteristic variation scale of the solutions is much 
smaller than that of the coefficients. Then, local solutions of 
the form 

i,fj ~ exp(ifc r r), (10) 
lead to a dispersion relation for the radial wavenumber k r , 



2 ) 



(11) 



By introducing the effec tive t ransverse wavenumber fe tr , de- 
fined by Bildsten et al. (1996) as 



k 2 



(12) 



the dispersion relation may be re-written in the more useful 
form 



, 2 2 2 _ f, 2 2 



2 ) Ov 2 -- 2 ), 



(13) 



The value of k r for given u and r, calculated using this ex- 
pression, determines the local character of waves at the ap- 
propriate frequency and location within the star. Inspection 
of equation (hoj) shows that real values (k 2 > 0) correspond 
to propagative regions, where the waves oscillate spatially, 
whilst imaginary values (k 2 < 0) correspond to evanescent 
regions, where the waves grow or decay exponentially in am- 
plitude. The k 2 = curves in the (r, u> 2 ) plane, defined by 
the roots of the right-hand side of the dispersion relation, 
separate these two types of region, and therefore correspond 
to the reflective boundaries discussed in the introduction. 
These boundaries, of fundamental importance when trap- 
ping zones are considered, are examined in the following 
section with the aid of propagation diagrams. Note that this 
k r — definition of the reflective boundaries formally vio- 
lates the assumption used previously to derive the solutions 
(p"o|); however, this violation will have little effect on the 
positions of the boundaries, and is not important at a qual- 
itative level. 

The remainder of this section is left to a discussion 
of the effective transverse wavenumber fc tr , since, as will 
be demonstrated subsequently, this quantity can be pivotal 
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Table 1. The physical parameters of the stellar model considered 
throughout. 



M/Mq 
7.00 



R/R Q 
3.16 



L/L Q 
1.76 x 10 3 



Teff/K 

21,000 



Z 
0.02 



in determining the trapping conditions at the stellar sur- 
face. In the case of plane waves in an infinite, plane-parallel, 
stratified medium, fc tr may be regarded as a free parameter; 
however, in the case of a spherical configuration, its is con- 
strained to assume values permitted by equation ( |l^ ) . These 
constraints arise from to transverse boundary conditions ap- 
plicable to waves propagating in horizontal (i.e., non-radial) 
directions; in a non-rotating star, they are equivalent to the 
requirement that solutions are invariant under the periodic 
transformations 6 — > 6 + 2n and cj> — > cj> + 2n, and lead to 
the familiar result (e.g., Unno et al. 1989) 



2 _ 1(1 + 1) 

fni — n — 



(i/ = 0). 



(14) 



When significant rotation is introduced, both of these 
requirements still hold, but an additional constraint in 9 is 
introduced as a consequence of the variation of the Coriolis 
force with latitude. Such variation means that, for v > 1, 
waves near the equator which are propagating in the lati- 
tudinal direction become evanescent when \n\ > 1/v; subse- 
quent reflections lead to the trappin g of these waves within 
the so-called 'equatorial waveguide' (Gill 1982). The result- 



ing horizontally-standing waves, whose angular dependence 
is described by the Hough functions Q] n (^; v), are oscillatory 
in latitude between the waveguide boundaries at /i = ±1/^, 
and evanescent elsewhere. With increasing v, these bound- 
aries converge towards the equator; for significant rotation, 
the constraints on fc tr therefore become dominated by the 
approximate requirement that an integer number of half- 
wavelengths in latitude fit between the waveguide bound- 
aries. This requirement is manifest in the a symptotic ex- 
pression for A; m found by Bildsten et al. (199f), which, when 
substituted into equation (|l2j), gives 

,2 (21, - 1) V 



(v » 1), 



(15) 



where is the number of latitudinal nodes exhibited by the 
appropriate Hough function between the waveguide bound- 
aries. This latter quantity is independent of v for prograde 
(m > 0) and zonal (m = 0) modes, whilst it increments 
by 2 as v is increased beyond unity for retrograde (m < 0) 
modes, due to the introduc tion of an addition a l pair of lat- 
itud inal nodes at v = 1 (Lee & Saio 1990; Lee & Saic 
199? |). f urthermore 



1, = I — \m 
sociated Legendre polynomials P{ 
over -1< n < 1, and 6r(M;0) = 



for v = 0, since the as- 
(fi) exhibit I — \m\ zeroes 



P/"(/i). These properties 
mean that in the above asymptotic expression (113) may 
be written in terms of I and m as 



'Mr 



(21 - 2\m\ ± 1)V 



(v » 1), 



(16) 



the plus sign being chosen for retrograde modes (m > 0), 
and the minus sign for prograde and zonal modes (m ^ 0). 
A comparison of this result with equation ([l4]) indicates that 
the permitted values of fctr for v 3> 1 can greatly exceed the 
corresponding ones in the non-rotating case, especially for 
small I ml. 



100 r 



10 



0.1 r 



0.01 h 




- 10 



6.0 5.5 
log T (K) 

Figure 1. The propagation diagram for I = 4, m = — 1 modes in 
the 7 Mq stellar model, plotted in the (log T, u) 2 ) plane to empha- 
size the outer regions of the model; the period 14 corresponding to 
the frequency is shown on the right-hand ordinate. Hatched areas 
correspond to propagative regions where waves have g-mode (/) 
or p-mode (\) character, whilst evanescent regions are blank. The 
dotted horizontal line shows the position of the trapping cut-off 
frequency wt- 



The notable exception to this discussion is the case of 
the prograde sectoral modes (m = — I), which in the limit 
v 2> 1 are transformed into equatorially-trapped Kelvin 
waves. Such Kelvin waves have a e xponential latitudinal de- 
pendence at small \x described by ( pill 198*2 ) 



exp 



r>2 2 



(17) 



indicating that they should be considered evanescent in the 
latitudinal direction even at the equator. Therefore, the con- 
straints on fct r are dominated by the periodic boundary con- 
dition in cj>- the exp(imei) azimuthal dependence of the so- 
lutions CPj — H) t nen gives the transverse wavenumber for pro- 
grade sectoral modes as 



(v > l,m = -I), 



(18) 



which can also be derived using the a symptotic expression 
for Xim found by Bildsten et al. (199£) for these modes. 



4 WAVE TRAPPING AND LEAKAGE 

As was demonstrated in the preceding section, the charac- 
ter of waves within a star is determined by the local radial 
wavenumber, so that positive and negative values of k 2 . can 
be identified with propagative and evanescent regions, re- 
spectively. An indispensable diagnostic tool for visualizing 
the location and extent of these regions, over a range of 
frequencies, is the propagation diagram introduced by Scu- 
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100 = 



10 



0.1 r 



0.01 h 




6.5 6.0 5.5 
log T (K) 

Figure 2. As for fig.[j], but rotation has been introduced at an 
angular frequency f2 = 8.04 X 10 rad s _1 . The azimuthal order 
m is -1. 



flaire ( 1974} ), in which the (r,Lo 2 ) plane is divided into zones 
over which the sign of k 2 is constant. 

Figure [l] shows the propagation diagram for I = 4 modes 
in a typical (non-rotating) early-type star; the logarithm of 
the temperature T has been adopted as the abscissa, rather 
than the radius, to emphasize the outer regions of the star. 
Regions in the (log T, uj 2 ) plane where waves are propagative 
(k 2 > 0) are hatched, whilst evanescent regions (k 2 < 0) 
are blank. Values for the Brunt- Vaisala, frequency N and 
adiabatic sound speed c a throughout the star, required for 



the evaluating k 2 using the dispersion relation 



aquin 



have 



been taken from a 7 Mq ZAMS stellar model, calculated by 



Loeffler ( 



199c 



by 

whose parameters summarized in table hi 
the model extends out to the photosphere at optical depth 
t = 2/3, where the temperature T = T e g = 21, 000 K cor- 
responds to an early-B spectral type. Equation (|l4|), which 
is appropriate in the non-rotating case, has been used to 
calculate fctr- 

In this figure, the type of hatching used to show prop- 
agative regions delineates between waves with p- and g-mode 
characters, the former occurring when both parenthetical 
terms on the right hand side of the dispersion relation ( |l3| ) 
are negative, and the latter when both terms are positive. 
The division of the diagram into relatively distinct p- and 
gr-mode propagation regions is characteristic of early-type 
stars. This division arises due to the fact that k 2 T c 2 diverges 
at the origin (due to the 1/r dependence of fctr) and is rela- 
tively small at the surface, while N 2 is approximately zero 
in the convective core (logT <; 7.3) and relatively large at 
the surface due to the steep stratification there. The promi- 
nent 'well' at log T w 4.6 indicates the presence of a thin 
convective region (N 2 < 0) due to Hen ionization, whilst 
the smaller well at logT ~ 5.3 is due to the metal opacity 
bump responsible for K-mechanism pulsation in early-type 
stars. 



As mentioned in the previous section, the k 2 = curves 
which separate propagative and evanescent regions corre- 
spond to the reflective boundaries required for wave trap- 
ping. Inspection of fig. |l| shows that, for gr-modes with fre- 
quencies greater than the trapping cut-off Ut, where uj 2 « 
1.9 x 10 -9 rad 2 s -2 is shown in the figure as a horizontal 
dotted line, there exists an extensive trapping zone formed 
by a pair of reflective boundaries, one at the edge of the 
convective core and the other in the envelope at lower tem- 
peratures. In contrast, for modes with frequencies below uj t , 
waves are propagative even at the surface of the star, and 
the outer reflective boundary required for the formation of 
a trapping zone does not exist. 

Strictly speaking, the term 'mode' is not appropriate in 
such circumstances, since stationary waves will not be estab- 
lished by repeated complete reflection. However, all waves 
at frequencies below the cut-off are evanescent in the con- 
vective region at logT ~ 4.6. This region, with a width of 
approximately 0.16 percent of the stellar radius, behaves like 
a partially-reflecting barrier to waves incident from the in- 
terior; some fraction of the waves will leak through the bar- 
rier and thence propagate unhindered to the surface, where 
they are lost from the star, whilst the remaining reflected 
fraction will contribute to the establishment of 'somewhat- 
stationary' waves interior to the barrier. Within the adia- 
batic approximation, these waves must decay exponentially 
in amplitude with time to compensate for the the energy 
lost through leakage, but will still exhibit a dis crete eigen- 
frequency spectrum. Shibahashi & Osaki (1976), when con- 
sidering a similar situation for high-frequency g-modes in 
evolved early-type stars, drew a useful analogy with virtual 
levels in the potential problem of quantum mechanics; there- 
fore, it seems appropriate to refer to such partially-trapped 
waves as virtual modes. Whether virtual modes can actually 
be self-excited in a star depends on the balance between the 
input of vibrational energy from a suitable driving mecha- 
nism, and the loss of vibrational energy associated with the 
leakage; non-adiabatic calculations are required to answer 
such a questions. 

The trapping cut-off frequency ui t , which separates the 
leaking virtual modes from the fully-trapped 'traditional' 
modes, is given by the smaller root of the dispersion relation 
( p^| ) at the stellar surface, namely 

-' 12 J ' ., (19) 



,2 2 1 

K tr c s 



where R is the stellar radius. This expression demonstrates 
the pivotal role of the effective transverse wavenumber kn, 
discussed at the end of the preceding section, in determining 
the trapping condition at the surface. In the non-rotating 
context, fctr can be eliminated from this expression through 
use of equation (jlij) to give 

Z(Z + l)c s 2 



2 

U t = 



(20) 



for n = 0. 

When the effects of rotation are included, the more gen- 
eral expression ( |i"2| ) for kti must be used in evaluating the 
sign of k 2 using the dispersion relation (|l^) . However, prop- 
agation diagrams may be constructed and interpreted in ex- 
actly the same manner as the non-rotating case. Figure [| 
shows the propagation diagram for the 7 Mq stellar model 
considered previously, but with rotation included at an an- 
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gular frequency Q = 8.04 x 10 -5 rad s~\ which is half of 
the critical rotation rate for the star, and corresponds to a 
period of 21.7 hours. The effects of the rotation on the equi- 
librium stellar structure having been neglected. Calculation 
of the eigenvalue A; m in equation ( |l2j ) , for each frequency or- 
dinate value in the (log T, uj 2 ) plane, was accomplished using 
Townsend's (1997) implementation of th e matrix eigenvalue 
method presented by Lee & Saio ( L990[ ). This method cor- 
responds to the spectral expansion of Hough functions in a 
truncated series of associated Legendre polynomials of the 
same azimuthal order m; 100 expansion terms were used 
throughout the calculations, a value deemed to provide suf- 
ficient accuracy since a similar calculation with 200 terms 
produced no numerical change in the results. An azimuthal 
order m = — 1 was adopted, so fig. ^| should be taken as 
appropriate for modes with (l,m) = (4, —1). 

Inspection of this figure shows that the trapping cut- 
off is significantly larger (wf ~ 8.0 x 10 -9 rad 2 s -2 ) than 
in the non-rotating case. This is a direct consequence of 
the influence of rotation on fct t ; at low frequencies where 
uj < 2Q, v > 1, and fc tr can assume large values, as discussed 
in the preceding section. The appropriate expression for u)t 
in rotating stars is given by 



Az m c s 



(21) 



although this should be regarded as formal, since it must 
be remembered that A; m is itself a function of uj through 
its dependence on the parameter v. However, in the limit 
fi 2 ^ c 2 /r 2 (at the surface), this expression will have solu- 
tions corresponding to v 3> 1, and thus the asymptotic ex- 
pressions ( [16|Jlq ) found previously may be used in the place 
of the general expression (jl9|) for fc tr . Solving the resulting 
equations for ujt then gives 



2^(22 - 2\m\ + l)c s /r 
2fi(2Z - 2\m\ - l)c s /r 



2 2/2 

i c B /r 



(m > 0) 

(-1 < m 0) 

(m = -I) 



(22) 



for Q 2 c 2 /r 2 . Applying the middle expression to the 
7M Q model for Q = 8.04 x 10~ 5 rads~\ and c s /r = 
9.76 x 10 -6 s _1 at the surface, leads to the asymptotic value 
uJt = 7.85 x 10" 9 rad 2 s -2 for (l,m) = (4,-1) modes, 
which is in reasonably good agreement with the value uj 2 ~ 
8.0 x 10" 9 rad 2 s~ 2 shown in fig. 0. 

The above expressions, when compared with equation 
(pp|), demonstrate that the effect of rotation is to increase the 
trapping cut-off uj t for all but the prograde sectoral modes; 
these latter modes will exhibit a smaller cut-off in rotating 
stars than in the non-rotating case, due to their transfor- 
mation into Kelvin waves discussed previously. This result 
is interesting in light of anecdotal observational evidence 
favouring prograde sectoral modes as the source of periodic 
line-profile variations in rapidly-rotating early-type stars. If 
such evidence can be substantiated a t a q uantitative level, 
as has been done by Howarth et al. ( 1998| ) for the rapidly- 
rotating pulsators HD 93521 and HD 64760, then it can be 
suggested that the bias towards prograde sectoral modes is 
due to the suppression of other types of mode, which will 
have a large values of uj t at rapid rotation rates and therefore 
preferentially leak from the star without self-excitation. 



5 EIGENFREQUENCIES 

In addition to its influence on the trapping cut-off frequency 
UJt, rotation modifies the eigenfrequencies and eigenfunctions 
of individual modes through its influence on the positions of 
trapping boundaries; this can be anticipated from the ap- 
pearance of A; m in the pulsation equations (^-Q). To evalu- 
ate the modified eigenfrequencies at a qualitative leve l, the 
asymptotic technique developed by Shibahashi (1979) and 



Tassoul (198C) may be adapted using the traditional ap- 
proximation to given expressions appropriat e for r otating 
stars. Using such an approach, Lee & Saio (|l987e|) found 
that low-frequency g-modes trapped between the boundary 
of the convective core (r = r c ) and the surface of a rotating 
early-type star have eigenfrequencies uj n given by 



(n + 7? e /2 - 1/6)tt 



\N\ 



dr. 



(23) 



where rj e is the effective polytropic index at the surface, and 
n is the radial order of the mode. This expression is not 
strictly appropriate in the current context, due to the fact 
that the outer reflecting boundary for trapped modes in figs. 
|l| and ^ occurs at r < R; furthermore, the presence of the 
convective region at logT m 4.6 has been neglected. How- 
ever, the form of the expression demonstrates that the fre- 
quencies of individual modes share the same A; m -dependence 
as the trapping cut-off uj t in equation ( pl[ 

It can therefore be suggested that a (/-mode which is 
trapped in a non-rotating star will remain trapped once ro- 
tation is introduced, since the effect of rotation is to scale 
both sides of the trapping condition uj n > uj t by an equal 
amount. Disregarding the possibility of avoided crossings 
(Aizenman, Smeyers & Weigert 1977; Lee & Saio 1989), a 
more general hypothesis may be put forward that, ceteris 
paribus, the set of radial orders {n} of the (/-modes which 
are trapped in a non-rotating star will remain invariant un- 
der the influence of rotation, with a similar result applying 
to the virtual modes. The hypothesis can be supported with 
an analogy drawn to atomic energy levels under the influ- 
ence of a magnetic field; even though the levels are distorted 
by the action of the Lorentz force (which, like the Coriolis 
force, can be expressed as a velocity cross-product), the set 
of discrete states which are bound is invariant under the ac- 
tion of the field. However, numerical calculations should be 
employed to test this hypothesis rigorously. 



6 DISCUSSION 

An important caveat regarding quantitative interpretation of 
the results presented previously is that atmospheric layers 
above the photosphere at r = 2/3 have been disregarded. 
This is justifiable if fc tr is constant throughout these layers; 
however, such a situation is unlikely to be realized, since 
even in the case of isothermal trans-photospheric regions 
where the sound speed c s is constant, fc tr ~ 1/r due to the 
spherical geometry. As a consequence, waves which are for- 
mally propagative at r — R may leak out to some radius 
r > R and thence be reflected back towards the interior, 
leading to complete wave trapping at frequencies below the 
cut-off UJ t . 

To address properly this issue of trans-photospheric re- 
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flection, it is necessary to relocate the nominal outer bound- 
ary of the star to a radius at which it is guaranteed that 
no reflected, inward-propagating waves will occur. In the 
context of the linear and adiabatic approximations adopted 
herein, this guarantee can only be made of the outer bound- 
ary is located at infinity. However, once non-adiabatic ef- 
fects are considered, it is possible that strong radiative or 
non-linear dissipation above the photosphere can lead to the 
effective absorption of all outward-propagating waves, with 
no reflection and subsequent trapping. Such a situation is 
analogous to the core-absorption of in wardly-propagating 
envelope p-modes found by Osaki (1977), and can be treated 
by placing the outer boundary at the base of the dissipative 
region (which may be close to the photosphere). 

Similar arguments concerning the absence of trans- 
photospheric reflection can be mad e for systems with stel- 
lar winds. Owocki & Rybicki (1986) found that, for a line- 
absorption driven wind, any wave-like disturbance which 
reaches the sonic point r 3 can never propagate back to 
smaller radii, due to the non- linear interaction between the 
wave and underlying mean flow. In this case, the outer 
boundary can be located at r s ; however, quantitative treat- 
ments are problematical, since the pulsation equations must 
be revised to take the sub-sonic wind regions into account. 

In spite of these difficulties, the results of this work are 
valid on a phenomenologi cal level, and may be of particular 
relevance to the 53 Per (smith & Karp 1976) and slowly- 
pulsating B (SPB; Waelkens 1991) classes of variable stars, 
which are unstable to low-frequency g-mode pulsation due 
to the metal opacity bump at logT w 5.3 (Dziembowski & 
Pamyatnykh 1993; Dziembowski, Moskalik & Pamyatnykh 
1993); the observed periods of these stars are typically 1 — 3 
days, which is the same order of magnitude as the trapping 
cut-off depicted in figs. ^ and []. As indicated previously, 
the global self-excitation of a given pulsation mode in one 
of these stars (or, indeed, any other type of star) depends 
on the competitive interplay between excitation and damp- 
ing mechanisms, which, respectively, pump energy into and 
remove energy from the pulsation at each point within the 
star; self-excitation will only occur if the net contributions 
from the former outweighs the net deductions from the lat- 
ter. For the purposes of the following discussion, both of 
these generic energy-transfer processes may be classified as 
either 

(a) non-adiabatic, where the transfer arises through per- 
turbations to the specific entropy, corresponding to the oper- 
ation of a Carnot heat engine which converts between ther- 
mal and mechanical (wave) energy within a given region, 



(b) advective, where the transfer arises through the non- 
zero divergence of the wave flux, corresponding to a net flow 
of mechanical energy through the boundaries of the region. 

The opacity mechanism operative in 53 Per and SPB stars 
is thus a non-adiabatic excitation mechanism |(a)[ whilst the 
energy loss associated with wave leakage at frequencies be- 
low the trapping cut-off uj t m ay be identified as an advec- 
tive damping mecha nism 



grangian pressure perturbation Sp tends to a limiting value 
at the surface; this corresponds to the ab initio restriction 
that all waves are evanescent at the surface, and thus com- 
pletely trapped within the star. Hence, any contributions to 
advective damping arising from wave leakage are neglected, 
which might lead to incorrect results for the overstability of 
modes at frequencies below io t . 

However, it must be stressed that the last point is some- 
what formal, if modes are stabilized by non-adiabatic damp- 
ing well before the frequency is low enough for leakage to oc- 
cur; whilst advective damping might enhance the stability 
of the virtual modes, it will have little influence in deter- 
mining which (trapped) modes are unstable in a star. In the 
case of the 53 Per and SPB stars, such a situation may arise 
due to the dominance of the opacity mechanism by radiative 
damping at lower frequencies. The latter is large for high- 
order (large-n) gr-modes, whose eigenfunctions exhibit many 
radial nodes in the stellar envelope; the sub-adiabatic tem- 
perature gradient in the radiative parts of the envelope will 
lead to significant thermal diffusion between neighbouring 
oscillating elements, which tends to suppress the pulsation 
(see, for instance, equation 26.13 in Unno et al. 1989 plus 
their accompanying text). These issues will be examined fur- 
ther in the next paper in this series. 

In a contrasting situation, where non-adiabatic damp- 
ing is less important at low frequencies, the overstability of 
virtual modes will be determined by the relative strengths 
of non-adiabatic excitation and advective damping. If the 
latter is dominant, then g-modes will exhibit an upper limit 
in their variability periods which corresponds to the trap- 
ping cut-off; in the more rapidly-rotating stars (see, e.g., 
Aerts et al. 1999), this upper limit will depend, amongst 
other things, on m and the degree of rotation. In contrast, 
if non-adiabatic excitation dominates, then no upper period 
limit will be observed, since virtual modes will be excited 
in addition to trapped modes. Estimates of the strength of 
advective damping can be obtained using, for i nstan ce, the 
asymptotic approach presented by Shibahashi_( [l979|) . How- 
ever, as with the local analysis used in section W, such an ap- 
proach is only valid when the characteristic variation scale of 
eigenfunctions is much smaller than that of the underlying 
star. This restriction means that Shibahashi's approach may 
lead to poor results for those virtual modes with frequencies 
close to Wt; therefore, an examination of the importance of 
leakage-originated advective damping is deferred to the fol- 
lowing paper, where the pulsation equations are solved glob- 
ally using a numerical approach which does not suffer from 
the restriction discussed. 

Whilst a proper treatment of trapping, even in cases 
without the trans-photospheric reflection described above, 
adds a certain level of complexity to theoretical studies, 
it does open the way for asteroseismological studies of the 
near-surface regions of early-type stars. For instance, if an 
upper period limit is observed as described, the inferred 
value of cot may be used, in tandem with equations (p(i[-^2[), 
to calculate a value for the acoustic timescale 



r/c B 



(24) 



(b) , In the stability calculations of 



Dziembowski et al. ( 1993), who use the approach described 



by Dziembowski ( 1977 ) , the assumption is made that the La- 



in the region where the onset of wave leakage occurs. Since, 
for an ideal gas, the adiabatic sound speed c s is a function 
of temperature T alone, this timescale then gives an inde- 
pendent estimate of the temperature in the outer layers of 
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the star. Conversely, observations of variability attributable 
to virtual modes can confirm the existence of sub-surface 
convective regions arising in ionization zones, predicted by 
evolutionary models of early-type stars. The degree of wave 
leakage associated with a virtual mode depends on the thick- 
ness of these regions (which form the partially-reflective bar- 
rier necessary for the existence of virtual modes); therefore, 
it might be possible to obtain estimates of the thickness 
through measurements of the leakage rate. 



7 CONCLUSIONS 

The prime conclusion to be drawn from the work presented 
herein is that the complete trapping of low-frequency g- 
modes beneath the surface of early-type stars is not guaran- 
teed. This is especially the case in rotating stars, where the 
trapping cut-off frequency u>t can be significantly increased 
by the action of the Coriolis force for all but the prograde 
sectoral modes. The fact that the latter are more effectively 
trapped in rapidly-rotating stars than other types of modes 
may explain anecdotal observational evidence which points 
to their favoured excitation. As a consequence of the depen- 
dence of g-mode eigenfrequencies on the rotation rate, the 
hypothesis has been put forward that the set of radial orders 
{n} of trapped <?-modes is invariant under the influence of 
rotation. 

Stability analyses which contain the ab initio assump- 
tion of complete wave reflection at the stellar surface might 
be in error at frequencies below the cut-off u t . More rig- 
orous calculations can include the possibility of wave leak- 
age, by adopting a more physically-realistic outer mechani- 
cal boundary condition. Such calculations will reveal to what 
extent advective damping associated with leaking virtual 
modes might suppress the self-excitation of these modes. 
These points may be of especial relevance to the 53 Per and 
SPB classes of variable stars. 
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